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We calculate exactly the velocity and diffusion constant of a microscopic stochastic model of 
N evolving particles which can be described by a noisy traveling wave equation with a noise of 
order N^ 1 ^ 2 . Our model can be viewed as the infinite range limit of a directed polymer in random 
medium with N sites in the transverse direction. Despite some peculiarities of the traveling wave 
equations in the absence of noise, our exact solution allows us to test the validity of a simple cutoff 
approximation and to show that, in the weak noise limit, the position of the front can be completely 
described by the effect of the noise on the first particle. 
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Traveling wave equations such as the Fisher-Kolmogo- 
rov-Petrovsky-Piscounov (F-KPP) equation Q, H H] 

d t u = d 2 x u + u-u 2 (1) 

describe moving fronts |4| in a large number of problems 
in biology, chemistry and genetics. In physics, it appears 
in non-equilibrium statistical mechanics and in the the- 
ory of disordered systems [a,lf|- In typical cases, ut{x) 
represents the concentration at position x and time t of 
some chemical species or of individuals carrying a certain 
gene. It is well established that equations of the F-KPP 
type have a one parameter family of traveling wave so- 
lutions with u{— 00) = 1 and u(oo) = parametrized by 
their velocities v, and that the decay of the initial con- 
dition determines the velocity of the front. For localized 
initial conditions, the velocity is the minimal velocity al- 
lowed fmin = 2 30,0. 

When deriving traveling wave equations such as (TO 
from a given microscopic stochastic models d H 0, 0, 
HI Hii ] , one usually gets a noisy version of JTJ : 



dtu = d x u + u — u 2 + \feg(u) t), 



(2) 



where the additional term is proportional to a func- 
tion g(u) of the concentration and to a Gaussian noise 
l;(x,t) white in time and possibly correlated in space. 
The amplitude e represents the ratio between the micro- 
scopic and the macroscopic scales. Internal fluctuations 
due to the finite number of interacting particles gives 
lla g(u) oc y/u(l - u) . Although it is now 
established 11 1 that the presence of the noise term is suf- 
ficient to select a single velocity v e , and that 
this selected velocity v e tends slowly j^, [lfl in the limit 
e — > to the minimal velocity Vm\r, allowed by Q, it is 
still a theoretical challenge H Efl El El El to predict 
how v mm — v e vanishes with e or how specific properties 
of the noisy equation such as the diffusion constant D € of 
the front position behave in the small noise limit. Other 
types of noise have also been considered; for instance, 
taking g(u) oc u in (J2J would represent the fluctuations 
of a control parameter due to some external noise [lfllEfl |. 



Here, we consider a microscopic model which can be 
viewed as the problem of directed polymers in a random 
medium |2 lj where we take the infinite range limit in 
the transverse direction. As explained below, this prob- 
lem leads to a noisy traveling wave equation and, for 
one specific choice of the disorder in the directed poly- 
mer problem, we can solve the microscopic dynamics and 
calculate exactly the velocity, diffusion constant and all 
the higher cumulants of the position of the front. This 
exact solution allows us to test several approximation 
schemes used recently to attack the general problem of 
noisy fronts. 

The microscopic model we consider here is defined as 
follows: at each time step t we have N particles on the 
real axis at positions x\(t), . . . ,xpf(t). Given these po- 
sitions at time t, the new positions at time t + 1 are 
obtained by 



Xi(t 



1) = max \xi(t) 

1<j<N J 



(*)] 



(3) 



where the Sij(t) are independent random numbers gen- 
erated according to a given distribution p(s). One can 
think of — Xi(t) as being the ground state energy of a 
directed polymer of length t ending at position i in the 
"transverse" direction. Then J^Jl describes the infinite 
range case in the transverse direction as the directed 
polymer may jump at any step t from any site j to any 
site i and gain a bond energy — Sjj(t). It is easy to show 
that the points remain grouped and we want to know 
the velocity w C xact and the diffusion -D C xact of this cloud 
of points (or, equivalently, of its center of mass). 

One can associate a traveling wave equation to (@ by 
considering the proportion it t (x) of particles on the right 
of x: 



u t {x) 



1 

N 



E « 

Ki<N 



(Xi(t) 



(4) 



where 9(z) is the Heaviside function. Clearly, ut(x) is a 
decreasing function with ut{— 00) = 1 and ut(oo) = 0, so 
that ut(x) has the shape of a front. 
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For a given front configuration u t (x) , the random vari- 
ables {xi(t +1)} are uncorrelated and the probability to 
be on the right of x is given by 

Prob(xi(t+l) >x) = 1- Yl M^-^W), (5) 

l<j<JV 

where n(s) is the probability that Sij(t) < s: 

/*(*) = f POO ^' (6) 



For a given «t(x), the probability JSJ) is simply the 
front (lit-i-i(x)) at time t+1 averaged over one time step. 
One could rewrite (JjJJ using only the front variable Ut{x): 



(u t+1 (x)) = 1 - exp 



, (7) 



but one needs to be careful when /i vanishes and use 
the following prescription: When /i(x — y) = 0, then 
w t (y)/i'(x — y)/fi(x — y) should be set to if u t (y) = 
and to +oo if u t (y) =f= 0. 

The random positions Xj(t + 1) can be generated by 
solving (ut+i(xi)) = Zi, where z\, . .. ,zn are N indepen- 
dent random numbers uniformly chosen between and 1 
and the fluctuating front ut+i(x) is finally given by 



u t +i(x) = e^(u t+1 (x)) 



Ki<N 



(8) 



From ©, one can calculate how the fluctuations of 
Ut+i{x) are correlated. Writing 



ut+i(x) = (ut+i(x)) + —;= Vt+i(x), (9) 
v iv 



one finds for x < y 

{ m+1 (x) Vt+1 (y)) = (l - (« t+ i(a:)>)<« t+ i(v)) (10) 



Writing T)(x) = t;(x)-\/ (?y 2 (x)), Eq. ljUJ) becomes very sim- 
ilar to with <?(u) = — u). So far, (jJUHJl and their 
consequence (j^l llUll are exact, for arbitrary N. From (JSJ) 
one can also calculate higher correlations of 774+1 (x) and 
show that, for large N, the 77 become Gaussian. For ex- 
ample one can check that, up to terms of order 1/N, 

(mmmm) = {mm) {mm) + {mm) {mm) + {mm) {mm) 

(We used the simplified notation t]i = ij t +i(xi).) We 
should notice that the Gaussian character of the 774+1 (x) 
is a property valid only for large N and in regions where 



Ar(l-(u t+ l(x))) >1 and N(u t+1 {x)) > 1. (11) 



From now on, we will limit our discussion to the case 
where p(s) is a Gumbel distribution: 



p(s) = exp (— s 



(12) 



In that case, the full analysis of the problem becomes easy 
and we can calculate exactly the statistical properties of 
the front in the large N limit. From © and gumbel, one 
has //(s)///(s) = e~ s and Q becomes 



(u f +i(x)) = 1 - exp [-B t e x ] , 
where B t is defined by 



B t — N I e x u t (x) dx = } 

Ki<N 



E 



Xi(t) 



This definition and (jl3|l imply that, given B t , 



(B 



t+l/ 



00, 



(13) 



(14) 



(15) 



which means that the distribution of the random variable 
Bt+i decays slowly. 

The main advantage of the Gumbel distribution i|12|) 
is that the maximum of several Gumbel variables is 
itself distributed according to a Gumbel distribution. 
Therefore, for fixed Xj(i) in (PJ), having the Sij(t) dis- 
tributed according to the Gumbel distribution l|12(l im- 
plies that Xj+i(f) is itself a Gumbel variable. 

The Gumbel distribution is simple because B t in i|14|) 
is the only information needed to construct the front 
Wt+i(x) at time t + l. If one defines the position X t of 
the front at time t by X t = In B t , then the displacements 



AX t = X t+ i -X t = InBt+i - hi Bt 
are uncorrelated random variables given by 
1 1 



AX t = In 



yi(t + i) 



VN(t + l) 



(16) 



(17) 



where the yi(t + 1) = B t exp[— Xi(t + 1)} are independent 
and, using (|5l 171 H3^l . distributed according to 



p(y) 



(18) 



As the AX t are independent, the cumulants of the posi- 
tion X t are simply t times those of AX t , which can be 
calculated from its generating function. The identity 



(e- 5AX <) 



1 



T(5) 



du 



.5-1 



<exp(- M e AX *)), (19) 



and the fact that the yi{t) are independently distributed 
according to i(T%|) , lead to the following exact expression: 



1 



m Jo 



du u 



5-1 



dy 



(20) 
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For large N, the integral over u is dominated by the 
neighborhood of u — where 



ydy~ l+ulnu+(2j E -l)u+0(u 2 \nu) , (21) 



where 7_e = — r'(l) = — / + °° du e~ u \nu is the Euler 
gamma constant. Any term of higher order in u would 
give a correction 1/N to the final result. Using i|21|l into 
(|2Up. one obtains, for large N, 



ln(e- 5AXt ) 
S 

~ L 



S(L + lnL) 

r'(i 



]nL + 1 - 2 7£; - 



•<0 



r(i + j) 



(22) 



where L = In A. This leads to the following exact ex- 
pressions for the velocity and of the diffusion constant of 
the front described by (|7 HTU|) when p(s) is given by (jT^ : 



Woxact = hm 



= L + lnL- 



-Doxact = hm 

t — >oc 



t 



TT 

3L 



InL 

2 



IE 



L 



(23) 



Expanding ill'L'l) in powers of 5 gives also all the higher 
cumulants of the position of the front. 

Armed with these exact results, one can test the qual- 
ity of various approximation schemes [ToL |22| . 

A first approximation is to neglect the noise and to 
write a deterministic traveling wave equation for the 
front. However, if one replaces (u t +i(x)) by u t+1 (x) in 
Q, one obtains for (|12ll a meaningless front equation: 
starting with uq(x) = 6{— .t), one finds ii2(x) = 1 for 
all x. Equation (jTJ is meaningful only in the presence of 
noise. 

Another way of removing the noise is to assume that for 
each position Xi(t+1) in ©, all the Xj(t) are uncorrelated 
random variables chosen independently for each i. This 
leads to the following deterministic equation of a front 
propagating into an unstable state: 



u t+ i(x) = 1 



dy u t (y)p(x - y) 



N 



(24) 



This equation is very similar to Q and its velocity 
can be obtained using the usual method^]: looking 
for solutions of the form Ut{x) — exp[— j(x — v^t)] 
when Ut{x) <C 1, one obtains a function ^(7). For 
initial conditions which decay fast enough, the veloc- 
ity is fmcanficid = mm-y u(7). When p(s) is the Gum- 
bel distribution, one gets 11(7) = [In N + lnT(l — 7)]/ 7 
and, for large N, the minimal velocity is Wmeanfioid = 
foxact + 1 + o(l/ L). One could try to improve this result 
by using the cutoff approximation [1(1 Il7i in pi)l. 
but as ^(7) depends on N, it is not clear that the for- 
mula Av = — 7T 2 7 2 f" (7)/(21n 2 N) derived in [ljj can be 



applied. Trying to apply it anyway, one obtains the ve- 
locity fmeanfieid — 7r 2 /2 which is not closer to Uexact than 

^mcanfield- 

The cutoff approximation can also be applied directly 
on the evolution equations l|13l 1141) by setting u t +i(x) — 
(ut+i{x)) whenever (u t +i(x)) > X/N and u t +i(x) = 
otherwise. One can write a closed expression for the evo- 
lution of B t : 



B t +i — N e x (l-e 



x 1 1 _ -B % e 



dx, (25) 



where the position A t +i of the cutoff is given by 
A 



N 



= (u t+1 {A t+1 )) = 1 - 



By eliminating A t +i, one gets 



B t+1 =B t N 



1 - e _u 



du. 



In N U 



(26) 



(27) 



Using u cu toff = ln(B t+ i/B t ), we get, for large L = In A 
Cutoff = L + InL + l ~ lE L + (j) . (28) 

Comparing l|28|l and exact, we see that w C utoff gives 
correctly the leading orders in L with a discrepancy 
(lnL)/L, which is only slightly above what the cutoff 
approximation may predict anyway, as there is no reason 
to choose any particular value of A, and which is much 
better than the mean- field velocity obtained from 124|l . 

So far, all the approximations replaced the noisy dy- 
namics by a deterministic equation. Thus, they gave no 
prediction for the diffusion constant. We will now exam- 
ine approximations in which the system remains noisy. 

A first possibility is to consider that the evolution is 
given by JH1 E3 Dill 114(1 where the noise term in is 
exactly Gaussian, even outside the validity range Ulljl. 
In this case, B t +\ is Gaussian but with (B t +i) = 00 (see 
((ToTl l. The noisy front equation © becomes meaningless 
after a single time step. So, one cannot ignore that the 
noise is not Gaussian near the rightmost particle. 

The next approximation one can try is to put a cutoff 
into the system as in (|26|) . and to take into account the 
effect of all the fluctuations on the left of this cutoff. In 
other words, at each time step, we determine A t +\ by 
iffify . we set u t +i(x) — for x > A t +\ and we use © 
for x < A t +i with a Gaussian noise n t +i(x) correlated 
as in (|10J) . This leads to a Gaussian B t +i characterized 
by (B t+ i) given by the right hand side of or Bt- 
evolution and 

{B 2 t+1 ) - (B t+1 f = 2N [ At+1 dy (1 - e~ B ^ V ) e y 



x I dx e x e- Bte ~ 



(29) 
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For large TV, one easily gets (B t +i) ~ NB t \nN and 
(B 2 +1 ) - (B t+1 ) 2 ~ 2TV 2 J B 2 /A, giving for the diffusion 
constant of the front position X t 



D 



(B 



'+1/ 



(B 



t+i/ 



■h- (30) 



Compared to the exact result l|23[l. this has the wrong 
L dependence and also depends on the precise value of 
the cutoff A. The velocity u cu toff — D/2 obtained in this 
case is also not a better approximation than the cutoff 
velocity l|28l) . We conclude that the effect of the noise at 
the left of the cutoff in the present model is too small to 
explain the value of the exact diffusion constant (|23|) and 
can therefore be neglected. 

A last approximation is to keep as only source of noise 
the stochastic position x max of the rightmost particle [22^: 
we choose £ max using the distribution of the exact front 
dynamics and, on the left of x max , we use for u t +i(x) 
the average value of the front given that the rightmost 
particle is on x n 



u t +x(x) 




In other words: 



if x > x u 

h + ^j^(u t+1 (x)) if x < x„ 



(31) 



with Prob(g) = e- q 6(q), (32) 



and from (|14|l . 



B 



t+i 



B t q 



1) 



1 - e 



du. 



q/N 



(33) 



Notice that, in that approximation, (B t+ i) — oo due to 
the contribution of small q or, equivalently, of large x max . 
For q small compared to TV, one obtains 



B, 



t+i 



B, 



= TV 



- + In TV + 1 - 7b - \nq + o(l) 

q 



(34) 



Using the definition of the displacement AX t and 
the identity (|19(l . one can compute from l|34l) the gen- 
erating function (exp(—5AX t )} . Up to order 1/lnTV, 
one finds nearly the same result as in the exact solution: 
one just needs to replace 1 — 2-f E in by 2 — 2j E . 

This means that, up to order 1/lnTV, the velocity of 
the front in this approximation is shifted by the small 
amount 1/ In TV, and all the other cumulants are the same 
as in the exact model. 

In this work, we have shown how the problem of di- 
rected polymers with N sites in 1 + oo dimension can be 
reduced to a noisy traveling wave equation |0 |SJ. For 
one special choice of the bond disorder 1|12J) . we could 
calculate exactly the velocity and diffusion constant Ij23(l 
of the front, and even all the cumulants J22J). The rea- 
son which makes this case soluble is that, at each time 



step, the only information one needs to keep about the 
past is a single variable B t given by 114fl . This is simi- 
lar to what was observed recently for shocks in exclusion 
processes H3, 13 where, in certain cases, one can de- 
couple the evolutions of the position of the shock and of 
its shape. Comparing several approximation schemes has 
shown that, in the present case, the cutoff approximation 
[lH Hil IrR gives a good estimate of the front veloc- 
ity, and that the full large N fluctuations of the front can 
be obtained by considering the effect of the noise on the 
rightmost particle. The predominance of this rightmost 
particle might be related to some difficulties noticed in a 
previously studied growth model [26j where only the first 
cumulants of the heights were considered. 

The front described by 10 is peculiar because TV ap- 
pears both in the noise term and the traveling wave equa- 
tion itself and because neglecting the noise in JJJ) leads to 
an ill defined traveling wave equation. In the mean field 
approximation, one obtains a F-KPP-like front equa- 
tion (|24|l which still depends on N and its velocity di- 
verges like In N . These peculiarities make the problem 
considered here rather different from usual traveling wave 
equations and does not allow us to use the present exact 
solution to check the validity of the ln~ 2 N shift of the 
velocity and of the In -3 N dependence of the diffusion 
constant which have been suggested by a number of nu- 
merical calculations [22I . The approximations (cutoff 
or noise limited to the rightmost particle) successfully 
tested here should however be helpful to describe more 
standard front equations. 

Of course, from both the points of view of the theory 
of disordered systems and of the theory of noisy traveling 
wave equations, it would be interesting to attack the case 
of a general distribution p(s). A starting point could be 
to try to make a perturbation theory where the Gum- 
bel distribution would be a zero-th order approximation. 
No need to say that obtaining the In -2 TV correction to 
the velocity would require a delicate resummation of this 
perturbation theory. 
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